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Abstract 

We treat a very crude model of an exploding star, in the weak field approximation of the Brans-Dicke 
theory, in a scenario that resembles some characteristics data of a Type la Supernova. The most no- 
ticeable feature, in the electromagnetic component, is the relationship between the absolute magnitude 
at maximum brightness of the star and the decline rate in one magnitude from that maximum. This 
characteristic has become one of the most accurate method to measure luminosity distances to objects 
at cosmological distances^^]. An interesting result is that the active mass associated with the scalar 
field is totally radiated to infinity, representing a mass loss in the ratio of the "tensor" component to the 
scalar component of 1 to (2ui + 3) (ui is the Brans-Dicke parameter), in agreement with a general result 
of Hawking 3; . Then, this model shows explicitly, in a dynamical case, the mechanism of radiation of 
scalar field, which is necessary to understand the Hawking result. 

PACS numbers(s): 04.30.Db, 04.50.+h,04.80.Cc 

1 Introduction 

Although the problem of constructing an exact solution in the Brans-Dicke (BD) theory, that may be 
considered a generalization of Vaidya solution^, has been addressed in the past[5], this solution has not been 
found yet. That generalization must reduce to Vaidya solution corresponding to the external gravitational 
field of a radially radiating star, in the limit of GR, and also, in the limit that the radiation flux vanishes, 
we must be able to obtain one of the empty space solution of Brans 

The interest of having a solution in the BD theory associated with electromagnetic radiation is to use it in 
the construction of a model of an exploding star that includes the ejection of a shell of mater, electromagnetic 
radiation in a similar way as in a model in GR[7||S1|5], and scalar field radiation. This may be in consonance 
with a star model in which a highly evolved core implodes to nuclear densities while at the same time the 
outer most layers of matter are blown at high speed. This idea can be realized in the weak field approximation 
of BD theory using the method of Barros and Romero^UI to construct a generalization of Vaidya solution 
that it is matched through a spherical singular hypersurface to a satisfactory interior solution. 

The formalism of a singular hypersurface in BD theory has been developed and used in the descrip- 
tion of bubbles dynamics in extended inflation models of the Universe jl XI I12| . among others applications. 
While Barrabes et al|13j treat a singular hypersurface in BD theory in the Einstein representation, another 
authors [TT1 I12| develop the formalism of a singular hypersurface in the Jordan-Fierz (J-F) representation. 
We shall follow the last approach in which it is usually considered that the metric corresponds to the physical 
one. 

In the weak field method of Barros and Romero if the metric field equation is considered in the Einstein 
representation, it does not depend on the scalar field. As a consequence, the dynamics of the singular 
hypersurface can be obtained in a similar manner as in the corresponding model of GR. Of course, the scalar 
field and the metric in the JF representation will depend on the matter model assumed for the ejected shell, 



while the radiation involved in the process will include two components, one corresponding to the scalar field 
and the other to the electromagnetic radiation, both tuned with the evolution of their sources. Thus, once 
we have determined the evolution of the sources, we know the interior and exterior solutions and the content 
and variation of the radiated energy at infinity in both components. 

In section 2 we present a brief summary of the results of the model of Barros and Romero^D]. Section 3 
is devoted to present the results of the theory of surface layers in BD theory, useful in the weak field limit. In 
section 4 we discuss a general model of an stellar explosion by constructing the metric and the scalar field both 
in the interior and exterior regions of the singular hypersurface. We end the section with a derivation of the 
frequency shift of spectral lines in the background metric and the luminosity of the exploding object as seen 
by a distant observer at rest. In section 5 we treat a particular model by matching two Vaidya solutions, with 
time dependent masses m~ and m + , respectively, through a shell of dustjEJ], to obtain the solution of Einstein 
equation, necessary to construct the BD solution in its linearized form, for the same matter distribution. 
The time evolution of the shell of dust is completely determined by the GR matching conditions for the 
metric alone. The resulting system of equations of motion consists of two ordinary differential equations for 
four unknowns. To complete the system, we obtain, first, an equation that shows how the total rest mass of 
the shell changes due to the balance between the total energy per unit time incident on the shell, J_, and 
the total energy per unit time emitted by the shell, J + , both at the shell surface. Then, according to the 
general characteristic of the scenario in which we are interested in, we may consider that J~ <C J + , and 
therefore, J~ it is neglected in the energy balance equation. This assumption is equivalent to consider vnT 
as constant. Another assumption is to take, essentially, J + proportional to the radius of the shell through a 
constant, x, that it is characteristic of the radiation production mechanism within the shell. This choice is 
similar to that one made by Hamity and Gleiser[7], on different grounds, that leads to a satisfactory model 
of a stellar explosion, with characteristic parameters such as maximum luminosity and time decay from that 
maximum, corresponding to a Type la supernova |2j. This is precisely the scenario in which we propose to 
compute the output of scalar energy from the ejected shell. For this particular example we compute the 
effective active mass of the system as a function of an external observer time. This effective mass presents 
two components: the "tensor" component, related to the radiation of electromagnetic energy, and the scalar 
component, related to the radiation of scalar field. We end the section with a numerical integration of the 
system of differential equations from a given set of initial conditions. We show a comparative set of results 
for different values of the separation constant x, the Brans-Dicke parameter uj, and the initial active mass of 
the shell. One interesting result is that the active mass, associated with the scalar field, is totally radiated 
to infinity; for the given system of initial conditions the process takes place in a short time, well before the 
star reaches its maximum luminosity if the initial active mass of the shell is of the order of or greater than 
1CP 2 M Q . This represents a mass loss in the ratio of the "tensor" component to the scalar component of 
1 to (2u> + 3), in agreement with a general result of Hawking j^J. Then, this model shows explicitly, in a 
dynamical case, the mechanism of radiation of scalar field, which is necessary to understand the Hawking 
result. In the last section we summarize the main results of our work. 



2 The weak field limit 

In the weak field approximation the solutions of BD equations are simply related to the solutions of GR for 
the same matter distribution. Following Barros and Romero|10| the corresponding field equations are: 

9ab{x) = (G Q c/)y 1 g ab (Go,x) , (1) 

9ab{x) = r] ab + hah , Kb < 1 , (2) 

<j>{x) = <f>o + e(» , (3) 

Gab = &TrGoT ab , (4) 
8ttT 

De = 2^3 ' (5) 

The "bar" on top of a symbol means that it is considered in the Einstein representation. For instance, 
9ab(Go, x) is the solution of eq.|@J). The same symbol without the bar is in the J-F representation which it is 
usually considered as the physical one. Thus, the physical metric is obtained from eq.Q). The box operator 
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□ corresponds to the wave equation for the flat metric; T a b is the (matter) energy momentum tensor. In 
first order T a b = T a b] T is the trace of T a i,; e(x) is a first order term in the energy density; cj>o is an arbitrary 
constant which satisfies |0o| ^ Go = 0q • Actually, in order that BD theory posses a Newtonian limit 
this constant must be related to the Newtonian gravitational constant G bv|14| 



V2w + 4 

Essentially, in the weak field approximation the metric calculated from BD equations is quasi-conformally 
related to the metric calculated from Einstein equations for the same matter configuration. The term 
quasi-conform means that in going from an Einstein solution g a b(G,x) to the corresponding BD solution 
g a b{x), apart from the conformal factor [1 — Go we must replace G by the new w-dependent "effective" 
gravitational constant Go given by ©. 



3 Singular hypersurface in BD theory 

The procedure sketched in the previous section requires to find first the metric g a b- To this end we have to 
specify the tensor T a b corresponding to the matter. We choose to have a singular hypersurface of matter, 
E, that separates the spacetime, M, into the external region M + , with metric g^ b , and the internal region 
M~, with metric g~ b : M = M + USU M~. The regions M + and M~ are non-empties; this a necessary 
requirement to show that the weak field condition (J2J guarantees the condition in M. 

Before we propose the solutions g^ b and in the regions we have to consider the matching conditions 
at E that such solutions have to satisfy. Of course, the matching conditions of the metrics g + and g~ are 
the same as in GR|15|: 

8ttG S = (K+ -K~)-g tr{K+ - K~) , (7) 

where K + and K~ are the corresponding extrinsic curvatures of E as "seen" from M + or M~ in terms of 
the metrics g^ b , respectively; K + and K~ are tensors in E; g is the induced metric on S by the metric g of 
the spacetime M. Finally, S is the surface energy momentum tensor in S. 

The extrinsic curvature K + can be expressed in terms of its components in a coordinate basis {e^} in E: 

K = ^)Kt b , K+ = - ( h c a h d b N c , d ) | E+ , K = 5 c a + N a N c , (8) 

where ^V| s+ is the unit normal vector to S, pointing from M" to M + , with its components expressed in 
AI + ; In eq. JHJ ";" indicates the covariant derivative using the Riemannian connection of M + associated 
with g~^ b . Similar expressions correspond to K^. 

On the other hand, eq.JSJ is in the J-F representation. The corresponding matching conditions on the 
scalar field <j> through E are: 

e| E+ = e| s - , (9) 
while its normal derivative has a discontinuity given bv|12|: 

Rirtr 9 

[JV(e)] E+ - [JV(e)] E - — , (10) 

3 + 2u> 



where trS is the trace of S, 



and a similar expression for E . 

4 The model 

4.1 General characteristics of the model 

To construct a simple model of an stellar explosion we choose E as the history of an spherical surface with 
metric 

dsl = dr 2 - R 2 {r)dn 2 . (12) 
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Also, we consider M + filled with a coherent unpolarized radial flow of electromagnetic radiation represented 

by 

T+ b = ( P k a k b )+ , (13) 
where fc+ is a null vector. Then, the exterior metric solution of eq.@ is represented by Vaidya's metrical]: 



with 



k +a = 5* , ^- = -Anr 2 p+ . (15) 
d\i 

The scalar field in the exterior region satisfies eq.JSJ) with T = 0. An outgoing wave solution is: 

e+(z) = M. , (16) 



where f(p) is an arbitrary function at the moment. 
We consider a general model corresponding to: 

• A time dependent spherical shell of matter with an energy momentum tensor, compatible with the 
symmetry of the problem, given by 

S = rjv®v+p(v®v — g) , the basis {e_y} = {d/dr, d/d6, d/dip} , (17) 

where v is the velocity of a fluid element in E, r/ the surface energy density and p the isotropic pressure 
within the shell. 

• The metric g~ b has the general form 

ds 2 = F 2 (£,r)df ~ H 2 (£,r)dr 2 -r 2 dil 2 , (18) 
where we have chosen the area variable r as one of the coordinates, for simplicity. 

• The metric g^ b in first order takes the form 

ds 2 + = U - 2g o^ + ^) _ Gq e +\ d f + (l - Go e+)(2dfxdr - r 2 dQ 2 ) . (19) 

• The interior of the shell is occupied by a central spherical body and, perhaps, electromagnetic radiation 
and traceless matter in a neighborhood of such that in that region, close to E, the scalar field 
solution of eq.© is 

r r 

i.e., it is the sum of a spherical wave moving towards decreasing r plus a spherical wave travelling 
towards increasing values of r. The functions g\ and g-i may be determined by imposing initial and 
boundary conditions in M~ and on E. 

From eqs.Q, © and IjlUI) we know that the scalar field and the dynamics of the spherical shell are closely 
related; that dynamics is contained in the solution of eq.@. From eq.|0 we obtain 

fHr)]=9i[ar) + R(r)]+g2[m-R(r)} ■ (21) 

To apply the matching condition IjlUI) we need to know the components N^- in the J-F representation of the 
normal vector N, with the requirement N^ 1 > 0. Those components may be obtained from the normalization 
conditions: 

vv=l , vN = , N-N = -l , (22) 
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where the components v a — (X, R,0,0) has to be computed on both sides of E, separately. In particular on 
E~ the time component X~ = d^(r)/dr; similarly, on E + is X + = d//(r)/dr. From IL'l't we obtain 

iV+ r = -B , 7V+^ = X+ , N~ r = -{F/H)X~ , N~ £ = -(H/F)R . (23) 

The definitions of the symbols used in H23[) are: 

1 /2 

B = + |B 2 + l-^(2m+-/)| , (24) 

X+ = (1+ G e+)(i? + B)- 1 . (25) 

The expressions for X~ , i/(£,r), and F(^,r) can only be known symbolically and up to the first order in 
our approximation. However, this knowledge is enough for the purpose of the present work. In particular 
we have 

X- = I - (1/2)[R 2 -poo(iZ,0] , (26) 

where </?oo(£> is the first order term in the expansion of the (00) component of the metric. A straightforward 
calculation in first order gives 

N(e)\^ = -(X- +R)^ + (X- -R) 9 4 + ^^X- , (27) 



R v ' R R 2 

X+f Bf 
R R? 



N(e)\ E+ = ~^ L + ^t , (28) 



the (') indicates the derivative of the function. Computing the derivative of expression 121|) with respect to 
r we obtain 

g' 1 (X-+R)+g' 2 (X--R)-f'X+ = . (29) 
Therefore, from eas. l(77|) . lf28|) and J2HJ) we have 



N(e)\ s+ - N(e)\^=(X-+R)^ + (B-X-)^ . (3.M 



Finally, taking into account that in first order (B — X ) ~ R 2 , g [(X + R) — dgi/dr, and ea. (|10() we have 

dgi _ AnR( V ~2p) R 2 f 



dr 3 + 2w 2R 



(31) 



Equation (|31|) is valid in the first order of approximation, including the terms in R 2 . This last equation, 
along with eq. J7J), complete the description of a general model. 

4.2 Frequency shift and luminosity 

Consider 1 in M + an observer at rest at spatial infinity and let t measure proper time along its world-line; the 
observer 4-velocity u = (d/dt) has components it M = 1; u? = with j = r,6,4>. Electromagnetic radiation 
emitted at the surface of the shell with characteristic frequency uj e = k ■ v(= X + ) will be received at spatial 
infinity with frequency uj r = k ■ u(= 1). Thus, 

1 + G 0£ + _ dt, (32) 



" r R 



R 2 + l-G ^4^- 



Va ~ dr 



We assume that (2m + — /) > 0. Notice that (u) e /ui r ) — > oo (infinite redshift) for R — ► Go (2m + — /) and 

it < o. 



1 From now on we choose units such that G = 1; i.e., Go = (2 + 3)/(2 u + 4). 
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The total electromagnetic energy radiated by the shell in all direction per unit time (the luminosity of 
the system) measured by an observer at infinity is given by 



L = lim 4tt r z (T ab u a N b ) =-m + /X 



,+ / Y+ 



(33) 



where N is the unit normal vector to the surface of constant r and u the 4-velocity of an observer at rest on 
that surface. The total electromagnetic energy £ em radiated in a time interval At is given by 



f — 



A/ 



Ldt = m + (t) - m+ (r + At) 



(34) 



In the case of a collapsing shell (R < 0), ea . I|33(l shows that L — > faster than m + as R — > Go (2m + — /) 
due to an extra redshift factor. 



5 A particular solution 

To present a particular simple solution of the general model we assume that: 

• The metric g~ is also represented, in part, by a Vaidya solution, corresponding to the exterior gravita- 
tional field of a radially radiating spherically symmetric central body (the core of the exploding star), 
given by 

ds I = (l - 2G ° m ' M ) dv 2 + 2dudr- r 2 dil 2 ; (35) 



i.e, we also consider that in the interior of M , close to S, we have a coherent unpolarized radial flow 
of electromagnetic radiation represented by 

T ab = {Pkah)~ , (36) 

where k~ is a null vector, with 

k- a =6 a r , ^- = -^r 2 p- . (37) 
dv 

• The hypersurface E is the history of a spherical shell of dust: 

S = riv®v . (38) 

The equation of motion for the surface layer, as given by the matching conditions eqs.Q and ijHJ) are 

R(A-B) = 4irG a r)R 2 =: m (39) 
1 ( h , G Q m+ G X+m+ \ 1 ( h Go m~ GoX^m_\ 

b[ r + -1p —)-l[ R + -lP —) = 47rGor '' (40) 

where 

X+ = (R + B)- 1 , (41) 

X- = {R + A)- 1 , (42) 

B . + («» +1 _i^) ,/J , (43) 

A - + . (44, 

We have introduced the total mass mo of the constituents if infinitely dispersed and at rest^S]. We assume 
that 

R(t) > 2 Go m+ > 2 G mT > . (45) 



G 



The momentum densities (or energy current densities) along the unit normal vector N to £ [E + (or E _ ) 
refers to £ when it is considered as part of M + (or M~)] measured by a local observer at rest on S (4- velocity 
v = d/dr) are given by 

(rJV»)± = : [ /9 (fc.iV)(fe-u)] ± =p ± (X ± ) 2 . (46) 

To write (|46|) we have used the orthogonality conditions N ■ v — 0, v -v — 1, N ■ N = —1, where all the scalar 
products are computed with the metric g. The total energy per time unit incident on (J~), or emitted by 
(J+), the shell is 

J± =: 4nR 2 (T Nv )± = - m k X ± . (47) 

Taking the r-derivative of l|39|) to use the resulting equation in H40|) . considering R ^ 0, and the definitions 
H4 1144( 1. we obtain^ 

too = J~ - J+ , (48) 

which is an energy-balance equation. It is easy to show that equation and the definitions l|4"3*)l and 

give 

j rh to - to m , _ 

A = 1 , B = , to = to — to . (49) 

too 2R too 2R 

If the functions to + (t) and m~(r) are given, (|39[) and l|4(J|l represent a system of second-order ordinary 
differential equations for the unknowns R(t) and t](t). 

5.1 The equation of motion as a first order system 

We are interested in a supernova explosion as the scenario to study, in a specific example, the main charac- 
teristics of scalar radiation in the BD theory. To this end we shall make some further assumptions to simplify 
the equation system. Let us consider first eq. H48fl . Once the gigantic explosion of the star takes place its 
core may collapse so rapidly that it forms a sort of extremely compact (degenerate) matter. This compact 
object, which may be a neutron star 2 or a black hole, is referred to as a compact supernova remnant. It 
may also be present a diffuse supernova remnant as a consequence of the shock wave and ejected material 
expanding from this explosion, and the interstellar material it sweeps up along the way. In our very crude 
model of the explosion, part of the interior of the shell is represented by a radiating compact object of mass 
to~(t) whose action on the rest mass of the expanding shell is represented by J~ in l|48|) . It seems then 
reasonable to assume 3 that J~ <C J + and neglect J~ in (|48|l . This assumption is equivalent to consider m~ 
as a constant. Thus, eq. (|48() becomes 

too = - J+(= - to+ X+) . (50) 

In consonant with the last assumption we assume that the rate of scalar radiation, that originates at E, 
to the interior and exterior of the shell of dust, are equal; i.e.: 



dgi 
dr 



= V 



(51) 



According to l|21|) what this assumption actually means is that g2 is constant; i.e., we have in part of A/ - , in 
first order, an incoming scalar wave that originates at E, plus a static scalar field associated with the central 
body. 

Let us introduce now the function a(r) by 

a(r) = to/toq . (52) 



2 At the moment of a neutron star's birth, the nucleons that compose it have energies characteristic of free fall, which is to 
say about lOOMeV per nucleon. That translates to 10 12 K or so. The star cools off very quickly, though, by neutrino emission, 
so that within a couple of seconds the temperature is below 10 11 K and falling fast. In this early stage of a neutron star's life 
neutrinos are produced copiously, and since if the neutrinos have energies less than about lOMeV they sail right through the 
neutron star and the surrounding matter without interacting, they act as a wonderful heat sink)l(i|. 

3 Since J~ and J+ are both functions of time, the least favorable instant is at r = 0. In the next subsection we shall see 
from the initial values of all the relevant variables that jf > 10 4 Jq. 



7 



In terms of the parameter a and using eqs. I(43|) . I|44|) . and (|49|l we obtain a first order equation for R(t) in 
the form|7j 

£22 1 . r m+ + m- G 2 rh 2 

R -a -1+Go ^ ^2 ■ ( 53 ) 

Similarly, from H5U[) and ((52(1 we have 

rha = arh(\ - aX + ) . (54) 

Replacing first B given in (|49|l into 1(41(1 and then the resulting expression for X + into (|54|l . we obtain 

rha(2aRR + 2a 2 R - G m ) = arh(2aRR - G m ) . (55) 

Eqs. 1)53(1 and 1(55(1 are two relations for the three unknowns, rh, a and R. Condition 1(45(1 now reads a > 0, 
2a 2 i? > G rh > 0. 

To find an explicit solution it is necessary to fix one of the unknown functions and then solve for the 
resulting differential equations for the other two. To guide our intuition let us replace the last obtained 
expression for X + into 1(54(1 and then use 1)55(1 to replace rh in the result. Thus, we obtain 

/+ 

(2aRR — G m ) . (56) 



2rhR 
We assume now 

J+ = 2 X rhaR . (57) 

The value of the constant x 1S correlated to the time scale of the process and it may be conjecture that it is 
a characteristic of the radiation production mechanism within the shell. For \ — we have a non-radiating 
system (a, rh — constants). The dimension of x equals the inverse of a square length; the adopted unit in this 
paper is Mg 2 . This choice is equivalent to that one made by Hamity and Gleiser[7], on different grounds, 
that leads to a satisfactory model of a stellar explosion, with characteristic parameters such as maximum 
luminosity and time decay from that maximum, corresponding to a Type la supernova. This is precisely the 
scenario in which we propose to compute the output of scalar energy from the ejected shell 4 . Finally, from 
(|Sl)jl and we obtain 

a = -x a {2aRR - G rh) , X = const. > . (58) 
Hence, we can write (|53f) and l|55ll in the form 

( i ~ rh + 2m- Glrh 2 \ 1/2 , . 

rh = - xm(2aRR + 2a 2 R- G Q rh) . (60) 

The system of first order differential equations, 1(581 1591 160(1 . that describes the motion of the shell, has 
to be considered in conjunction with the first order differential equation for the amplitude of the outgoing 
scalar wave, which according to ea. 1(31(1 . and assumptions 1(38(1 and 1(51(1 . becomes 

; . iTiR 2 ri R 2 f 

J = 9i = - 



i?(3 + 2uj) 2R 

rh R 2 f 
aR(3 + 2uj) ~ ~2R 



(61) 



The initial condition for /(r) may be obtained from the requirement that at the time of the explosion, 
r = 0, the exterior scalar field matches continuously to a BD static solution in the weak field approximation 
generated by a central body of active mass m^JOl; i- e -i 

t'-^U- < 62 > 

4 If we assume that the total electromagnetic energy per unit time, radiated by the shell, is proportional to its area times the 
absolute temperature T at its surface to the fourth power, we have that T 4 ~ rha/R, which is a decreasing function of time 
during the expansion period. 



8 



where fi and mf are the initial values of /(r) and m + (r) respectively. 

In order to compare the results of our model, in the electromagnetic mode, with observational data for a 
stellar explosion, such as a Type la Supernova|17|. for instance, we need to compute the light curve, L(t), 
from eq. 1|33[) . for values of < r < Ti, where n corresponds to the time when the maximum luminosity has 
decline in one magnitude. 5 From eqs. (|24ll . (|25|l . (|33[1 and 1|59|) the expression for the light curve becomes 



L(r) 



- rn 



R 



G m G / 



4 a 2 R 2 



R 



R 



1/2' 



Ggf 
R 



(63) 



It is apparent from ea. l|19|) that the Keplerian mass M (the active mass), measured by an orbiting object 
around the exploding star, in the present approximation is 



M = rn 



//2 



(64) 



i.e., the total active mass is decomposed into the sum of a "tensor" component m + and a scalar component 
//2^2j. Therefore, from ea. l|19fl the total scalar energy radiated in the time interval At is given by 



A^ = [/(r)-/(r + Ar)]/2 => £ = --/ 



(65) 



Thus, from eqs. (|61|l and (|65|l the total active mass radiated in the time interval [0, n], corresponding to 
the observer time interval [0,<i], associated with the scalar field is 



£<t>\ — 



rh R 2 f 
2aR{3 + 2uj) + ~AR 



(It 



(66) 



5.2 Numerical results 

We have performed a numerical integration of the system (|32I34I58I59I60I61I65|) , corresponding to the func- 
tions t(r), £ em (r), a(r), R{t), t%{t), /(t), and £^{t), with initial conditions, at r = 0: 

^ = 0, £ eml = 0, a\ - 1 = 1(T 3 , R t = 2 x l(T 2 i? = 9469.4M Q , 

rhi = (1(T 3 , 10~ 2 , 10- 1 ) M , fi = 2m+/(2w + 3), 8# = 0. 

The constant mT = 1M@, and \ was chosen to have three different values, corresponding to different 
cases. The initial value of the total energy per unit time emitted by the shell verifies 6 J 4 + > 3.5 x 10 37 ergs/s. 
Similarly, for these initial values and u> = 500, the initial velocity of the shell becomes Ri « 10000 km/ s. The 
value of ui = 500 it is generally accepted^] as the lowest bound implies from Solar system tests, although 
some more recent estimates I2Qj have raised this limit to u> > 40000, obtained using signal timing from the 
Cassini spacecraft [21] ■ The same authors 20 point out that this bound may be weaker on cosmological scales 
than in the solar system[22]- On the other hand, in a recent publication 23 it is shown that the parameter 
ui "runs" with the scale factor in a Friedman-Robertson- Walker metric in order that a BD theory serves as 
a successful model for dark matter - dark energy. In this model the value of the parameter oj is less than 
(—3/2) in the matter to scalar transition period; equal to (—3/2) in the BD scalar field dominated era; and 
greater than (3/2) in the scalar to acceleration transition period. In any case, in our BD model of an stellar 
explosion in the weak field approximation, the main characteristics of the model are quite insensitive to any 
value of oj greater than 2 since it results Go ~ 1 , except for the radiation rate of scalar field and the total 
energy radiated in the scalar mode as it can be seen from eqs. ijBl^jl. ljo5|) . and lIHBl) . 

For the given set of initial conditions we have computed the light curve L(t), for < t < ti, and the 
total radiated active masses £ em i and in the same time interval. In Table \T\we show a comparative set 
of results for different values of the constants x an d oj, and the initial values of m^; Ati =t m — t\, where t m 

5 Actually, the maximum luminosity should correspond to the B-bandUl but in our very crude model of a stellar explosion 
we do not have the possibility to specify the B-band from the model of radial electromagnetic radiation corresponding to the 
Vaidya solution. 

6 The corresponding value for a star like the Sun is J+ ~ 10 33 ergs/s. 
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X 


U! 


rhi 




At! 


ft 


^eml 


£<t>\ 


[M ]- 2 




[M@] 


[crgs/s] 


days 




[ergs] 


[orgs] 


5 x 10~ 24 


500 


10" 3 


3.53 x 10 43 


22.1 


15 x 10" 4 


1.04 x 10 50 


4.46 x 10 50 


10 -23 


500 


10" 3 


4.98 x 10 43 


15.8 


15 x 10~ 4 


1.04 x 10 50 


4.37 x 10 50 


25 x 10~ 24 


500 


10" 3 


7.89 x 10 43 


9.70 


15 x 10~ 4 


1.04 x 10 50 


4.22 x 10 50 


10 ~23 


5000 


10" 3 


4.98 x 10 43 


15.7 


15 x 10~ 5 


1.04 x 10 5() 


4.38 x 10 49 


10 -23 


50 


10" 3 


4.98 x 10 43 


15.7 


14.7 x 10" 3 


1.04 x 10 60 


4.26 x 10 51 


10 ~23 


5 


10" 3 


4.98 x 10 43 


15.6 


12.0 x 10" 2 


1.04 x 10 50 


3.04 x 10 52 


10 -23 


500 


10~ 2 


4.99 x 10 44 


15.8 


(i) 


1.04 x 10 51 


(i) 


10 -23 


5 x 10 4 


10" 2 


4.99 x 10 44 


15.8 


(i) 


1.04 x 10 51 


(i) 


5 x 10" 24 


500 


lO" 2 


3.53 x 10 44 


21.9 


(i) 


1.04 x 10 51 


(i) 


10 ~23 


500 


io- 1 


4.99 x 10 45 


15.6 


(i) 


1.04 x 10 52 


(i) 



'The active mass associated with the scalar field, /i/2, is totally radiated to infinity after a time 
tj <t\, well before the shell reaches its maximum brightness. For the values of tf, see the text. 



Table 1: A comparative set of results for different values of the constants x an d w, and m,; the other 
parameters and m~ = 1M , are the same in all cases. 



and t\ are the observer times at the maximum luminosity, L m , and when it has decline in one magnitude, 
L%, respectively 7 . By fi, we indicate the final value of / at observer time t\. When m, = 1O _2 M0, /(r) 
becomes equal to zero at a time tf « 716s; while if m, = 10 Mq, it results tf rj 1.5s; in both cases tf <C to, 
where to is the observer time at which R = 0; in all the cases t m < to- The reason for this fast approaching 
to zero of /(r) is that / is essentially proportional to rh, through the first term on the right hand side of eq. 
(16 11 1 . Finally, it may be important to mention that although we have worked with the exact Vaidya solution, 
for simplicity, the same numerical results are obtained if we drop in the equations all the terms quadratic in 
rh. 

6 Final comments 

We have treated a very crude model of an exploding star in the weak field limit of BD theory. However, the 
model presents aspects in the electromagnetic energy radiated to infinity that resembles some characteristics 
data of a Type la Supernova. The most noticeable feature is shown in the first three lines of Table i.e., 
the value of the constant \ relates the absolute magnitude at maximum brightness and the decline rate in 
one magnitude from that maximum. This characteristic has become one of the most accurate method to 
measure luminosity distances to objects at cosmological distances We also notice from Table 1, that 
we may scale the value of L m by changing just the initial value m,, without affecting the decline rate At\. 
These values seem to be independent of the value of lo within a wide range. Finally, the total electromagnetic 
energy, radiated in the time interval [0, ii], is independent of the value of \ and u>; the total energy £$x, in 
those cases in which it is not limited by the initial value fi, it is larger the smaller the value of u>. 

On the other hand, the total active mass associated with the scalar field, as given by fi/2, is totally 
radiated to infinity. This process takes place in a time lapse considerably smaller than the time in which 
the star reaches its maximum brightness after explosion, if the mass of the shell is of the order of or greater 
than 1O -2 M0. This represents a mass loss in the ratio of the "tensor" component to the scalar component 
of 1 to (2cj + 3), in agreement with a general result of Hawking 3 . 
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